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Calculating  Thermal  Boundary  Layer  in  a  Flow  at 
Canfcreasible  Qaa 

toy 

L.  Zysina-Kol ozhen 

Discussed  la  an  approximate  se.dempirical  method,  which  allows  with  sufficient 
accuracy  to  calculate  a  lan inary,  transient  and  turbulent  zones  of  thermal  boundary 
layer  when  a  flow  of  compressible  gas  is  directed  around  a  surface* 

The  solution  of  many  technical  problems  is  connected  with  the  necessity  of  cal¬ 
culating  heat  exchange  of  a  surface  with  a  stream  of  compressible  gas  flowing  around 
it.  The  boundary  layer  originating  thereat  in  dependence  upon  the  flow  condition  and 
the  nature  of  velocity  distribution  along  the  surface  can  be  over  the  entire  span  of  the 
surface  either  la  mi  nary,  transient  or  turbulent,  or  on  sections  of  surface  -laminary, 
then  transient  and  finally  turbulent. 

In  this  report  is  introduced  and  approximate  semi  empirical  method,  allowing  with 
satisfactory  accuracy  to  calculate  all  three  zones  of  thermal  boundary  layer, origi¬ 
nating  when  a  compressible  gas  is  flown  around  a  surface.  The  nature  of  the  method 
lies  in  the  following. 

Ve  will  investigate  a  plane  flow  of  canpressible  gas.  It  i^known,  that  for  this 
case  the  Integral  ratio  of  energy  in  Dorodnitsyn  variables  acquires  the  farm 
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(sign*  corresponds  to  retardation  paramatere). 


It  can  be  written 
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where 
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Using  formula  (3)  it  is  possible  to  change  equation  (l)  into  farm  at 
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We  will  introduce  parameters! 

It  —  -77-  °r  ^r  * 
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and  will  assunm,  that  they,  changing  along  the  surface  around  which  flew  is  directed. 


do  synonymously  determine  all  characteristics  of  thermal  boundary  layer*  We  will  assume, 

that  in  expression  (8)  the  influence  of  the  longitudinal  temperature  gradient  in 

IFn&r-  _d&_ 

characterized  by  the  complex  UQ  ,  and  G^,  does  not  depend  upon  dx  and  determines 
only  the  effect  of  the  Re -number.  In  this  ease  the  value  Gj  will  be  identical  for  the 
flow  around  the  profile  and  for  the  flow  around  a  plate  as  well.  We  will  determine 
the  function  G^.  in  accordance  with  data  concerning  heet  exchange  of  a  plate. 

By  examining  the  experimental  date  (fig.l)  it  be  canes  evident  that  when  the 
physical  constants  are  referred  to  retardation  parameters  •  the  formulae  for  calcula¬ 
ting  heat  exchange  retain  the  very  same  form  as  for  the  case  of  a  noncoapresaible 
flow.  For  the  transient  zona  of  the  boundary  layer  the  effect  of  compressibility  la 
reflected  on  the  coordinates  of  the  beginning  and  ending  of  the  transition  (kh-^g, 

kh  and  the  development  of  the  transition  process  after  its  origination  take* 
end  * 

place  in  such  a  way.  that  lines  N(Rex)  remain  parallel  to  each  other  over  the 

entire  investigated  range  of  change  in  M  number. 

It  is  evident,  that  the  curves  can  be  approximated  by  a  family 

•Footnote  —  ■  -  During  the  calculating  all  .physical  constant  are  referred  to 
temperature  of  retardation. 
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In  which  the  value  a  B  and  n  are  different,  but  unchanged  for  each  flow  condition  in 
the  boundary  layer.  For  the  transient  zone  the  coeffieimt  B  appears,  as  shewn  by  expert' 
meats,  to  be  a  variable  value  .changing  with  the  change  in  value  of  Re  number,  corres¬ 
ponding  to  the  point  of  beginning  of  transition  Roy-y,K«>^»  ^Ms  value  is  determined  by 

the  initial  turbulence  of  the  flew  ^4  ;.  magnitude  of  the  temperature  factor  and  M* 

* 

number,  i.e.tht  value  B  retains  constancy  within  limits  of  each  concrete  experiment, 
but  can  change  with  the  change  in  initial  and  practical  conditions  of  the  process. 
Substituting  expression  (10)  in  the  integral  ratio  of  energy  for  plate 
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and  using  formula  (9)  it  is  possible  to  obtain  for  function  Gj  the  following  expression* 
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Fig.  l.Pe pend ence  Nbjf*  NCBe*)  at  various  values  of  M-nunber* 

1-0 .24 |  2-  0J*4|  3-C.52*  4-0.73*  5-  l»°7l  6-  1.43*  I-Nux=*0.297  Re^lazninary  bounda- 
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ry  layer  |  II-  Ita*-  0.0255  -  turbulent |  ^  1. 

We  will  determine  function  Z  in  such  a  way  that  for  a  plate  it  will  be  poaaible 
to  accept 
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Then  for  Op  is  obtained  a  formula. 


gt  a^'t' \  (14) 

which  is  perfectly  analogous  to  formula  for  an  incompressible  flow'll. 

*■*  • 

Introducing  formulas  (8).  (9)  and  (14)  into  equation  (7)  it  is  possible  with 
the  aid  ef  none  amp  lex  calculations  to  bring  same  into  the  following  fonu 

ft -  (15) 
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For  an  incompressible  flow  was  shewn  experimentally [2 .that  the  function.set 
in  square  parentheses 


I  t  ( m  -  ])/  2 It 

can  be  well  expressed  by  fornula 


(IG) 
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The  value  a  was  found  to  be  independent  from  the  temperature  gradient  and  is  deter¬ 


mined  only  by  the  condition  ef  flow  in  boundary  layer. 

We  will  assume,  tnat  for  a  compressible  flow  is  also  possible  to  write  an  analogous 
fornula.  If  a  does  not  depend  upon  the  longitudinal  pressure  gradient,  then  it.  appar¬ 
ently.  should  have  an  identical  Talus  as  during  a  flow  around  a  prof tie. and  during 
the  flow  around  a  plate.  Then,  by  using  fornula  (13)  it  ia  possible  to  determine  a  by 
formula 

r  t: 

«“(«  + 1)  [-^(l->i)*-!J  .  08)  • 


where  a  -  value  for  an  perturbed  flow. 

If  we  would  consider  flows,  characterized  by  condition  T^/T*  «  & »  const,  then 
it  is  apparent  that  for  each  flow  condition  a  will  be  a  constant  value  depending  only 
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upon  the  flow  condition. 


In  this  ease  equation  (15)  can  be  easily  integrated  and  allows  to  determine  <£** 
for  the  thermal  boundary  layeri 

»  m  X 

f  Q  \  »*rT  (v*)'"'1  (  r  a  \ 

V=  (t)  ~L~  I  J  W  +  'l  Wr'T*'  ~f.  09) 

Here  A^  corresponds  to  the  value  of  the  coefficient  A  in  formula  (li|)  for  land  nary 

boundary  layer,  if  the  calculation  is  made  for  the  transient  zone,  and  corresponds 

to  the  value  A  for  transient  zone  if  the  calculation  is  made  for  the  turbulent  part 

n  A ** 

O  T* 

of  the  boundary  layer*  Consequently  Re5£  *  — Zj~* —  *8  determined  for  coordinate®  of 
the  beginning  of  change  over  from  calculation  of  la  ml  nary  section  in  first  oas*  and 
coordinates  of  the  end  of  transition  fr cm  calculation  of  transient  zone  in  second. 

For  convenience  of  and  to  speed  calculation  value  a  far  each  condition  of  flow  in 
the  boundary  layer  can  be  calculated  by  formula  (18)  in  the  necessary  range  of  para¬ 
meter  changes  and  presented  in  form  of  curves  (fig.2). 

Using  formulas  (19),  (9).  (l6)-(l8)  it  is  possible  to  obtain  an  expression  for 
th.  calculation  of  local  values  of  heat  exchange  coefficient 

n“-  mf*  ^ 

XM 

Comparing  (20)  with  the  corresponding  formula  for  incompressible  flow^*1*  4* 

possible  to  detect,  that  when  all  physical  constants  are  referred  to  flow  temperature 

retardation  T*  there  is  s  formal  analogy  between  the  type  of  forrnla  for  calculating 
o 

intensity  of  heat  exchange  in  a  compressible  gss  flew  and  in  an  incompressible 
flow  Nu'z.  By  compering  these  formulas  it  is  easy  to  obtain  a  relation 
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It  is  evident  from  thi3  formula,  that  for  the  case  'of  a  flow  directed  around  a 
plate,  when  the  rate  on  the  outer  linlt  of  the  boundary  layer  is  equal  to  the  speed 
of  the  oncoming  flew  and  or  a  ^..respectively,  at  above  described  repre- 

snatstion  of  formulas  should  case  to  a  ratio 


As  already  mentioned,  this  formal  ratio  correspond*  well  with  experimental  date 


See  attached  page  for  Figure  2  (fa) 

Fig. 2. dependence  a=*(/)  far  1 ami nary  section 
of  boundary  layer  (a)*  for  transient  section 
(b)  end  for  turbulent  (o)  at  -various  1^/T*  i 

1-0,5 I  2-0.6*  3-0.7*  4-0.8*  5-0.9,  6-1.0* 

7—1*2*  8—1.4*  9—1*6,  10—1.8*  11—2.0 
When  calculating  the  coordinate  of  the  point 
of  the  beginning  of  change  over  kh^  can  be 
determftflfr  by  the  Dorodnitsyzv-Loyteyanckiy 
method roT  an  approximate  determination 
of  the  length  of  the  transient  zone  we  could 
recoumend  tha  following  simple  considerations. 
The  processing  of  experimovte  (results)  listed 
in  fig.l,  and  analogous  results  of  other 
authors  showed,  that  the  parameter  rx, char¬ 
acterizing  the  relationship  between  coordi¬ 
nates  of  end  Khead  and  beginning  of  transition  kje^g,  doas  not  depend  upon  M  number. 
This  property  of  the  parameter  rx  allows  to  determine  its  value  by  empirical  curves, 
obtained  for  an  incompressible  flow  and  mentioned  in  report |Vj. 

Designations 

-  density, corresponding  to  retardation  parameters*  U0,  TQ-  rate  and  temper¬ 
ature  outside  of  the  boundary  layer*  %  -temperature  of  wall. 
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